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Abstract 

We consider a wave packet of a charged particle passing through a 
cavity filled with photons at temperature T and investigate its localization 
and interference properties. It is shown that the wave packet becomes 
localized and the interference disappears with an exponential speed after 
a sufficiently long path through the cavity. 
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1 Introduction 

It has been known for a long time that thermal photons can substantially disturb 
electron beams in an accelerator It also has been suggested that the 

problem of the classical limit of quantum mechanics (at least the problem of 
decoherence) can be solved in the presence of a reservoir of photons . Some 
models of an interaction of a quantum system with a reservoir confirming the 
decoherent behaviour have been discussed in refs.|QQ[^. These models are 
treated as an approximation to the quantum electrodynamics (QED). 

We discussed QED and some approximations to QED in ref. Q . We pointed 
out a difficulty (resulting from the ultraviolet singularity) with some approxima- 
tions to QED when applied in a discussion of the decoherence. The singular part 
of the QED propagator does not contribute to the decoherence in formal semi- 
classical calculations. In this paper we discuss the regular part. We calculate 
the time evolution of the density matrix of wave packets in an approximation 
of a small charge e treating the particle dynamics in a semiclassical approxi- 
mation. We neglect the vacuum fluctuations considering the effect of thermal 
photons on the evolution of the density matrix. We show a decoherence effect 
of the thermal photons. The vacuum fluctuations lead to divergencies at short 
distances. After a renormalization the contribution of vacuum fluctuations to 
decoherence is negligible for a large time. 

The effect of photons on decoherence has been discussed first by Joos and Zeh 
[^j. The general arrangement for interference experiments in an environment 
has been considered in L. Ford discussed a change of the vacuum 
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fluctuations with a change of the environment's geometry (the Casimir effect) 
and its role in the interference experiments. In such a case the relative strength 
of the vacuum fluctuations rather than its intrinsic value (which is infinite) is 
relevant. The decoherent effect of the black body radiation has been studied by 
Stapp but we disagree with his results and approximation methods. 

We are interested in quantum electrodynamics at finite temperature T de- 
termined by the density matrix (the Gibbs state) 

p = [Tr (cxp {-I3Hb)))-^ cxp(-/3i7H) 

where = KT, K is the Boltzmann constant and H^^ is the Hamiltonian for 
the quantum free electromagnetic field. 

We compute the correlation functions explicitely 

G^(0,x;t,x')j7 =Tr(Aj(0,x)A,(t,x')p) = 

^/dk|k|-icos((x-x')k)<5*[(k) ^ (1) 
(cos(c|k|i)(i + (^exp(/3?ic|k|) - l) ^ - ^sin{c\k\t)^ 

where 

The term i?ic|k| corresponds to the zero point energy (of vacuum fluctua- 
tions) whereas ?ic|k|(exp(?ic|k|/3) — 1)^^ is the average energy of thermal photons 
with the wave number k. The vacuum fluctuation (noise) is a measurable ef- 
fect and in general cannot be neglected. The virtual photons corresponding to 
the vacuum fluctuations are not directly observable. Then, the thermal pho- 
tons are described by the Green's function Gth = — Goo (note that Gth is 
real whereas Gf3 and Goo are complex; in quantum field theory the imaginary 
part of the Green's function is related to the pair creation and annihilation, so 
subtracting Goo means that the processes of pair creation and annihilation are 
neglected) 

Gth{^, x', t),i = ^J dk|k|-i,5*[(k) cos ((x - x') k) cos (ct|k|) (exp il3hc\k\) - 1)"^ 

(2) 

Let us note that Gth determines the real Gaussian random field. Subsequent 
computations can be performed either in the Fock space or by means of the 
functional integration . We do not explain the equivalence of both methods 
here but refer to our earlier paper 0. 
For a small |x — x' 

Gth{x,^' ,t)ji ~ Gth(0,0, t)ji = Sji^hcn^'^ dkkcos{ckt) ^exp{Pfick) ~ 1^ 

(3) 

is approximately x- independent. 
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2 The semiclassical approximation 

We approach the semiclassical limit of the wave function in the standard way 
treating the electromagnetic field A as a classical field . Then, the quantum 
electromagnetic field is realized as a random field with the covariance (2) . A 
solution of the Schrodinger equation 

ihdti>{t, x) = ^ i-ihV + ^Atfi,{t, x) (4) 

with the initial condition tp = exp{jrW)(j) can be related to the solution of the 
Hamilton-Jacobi equation Wt 

dtWt + :^{VWt + -Atf = (5) 
2m c 

with the initial condition Wt=o(x) = W{x). We express ipt in the form 

V't = XfPt = exp{-Wt)(f)t 

Then, ipt is the solution of eq.(4) if and only if (pt is the solution of the equation 

if) 1 p 1 p 

dtcPt = ^Act>t - -{VWt + -At)V4>t - ^{AWt + -divAt)cl>t (6) 

with the initial condition 4>. In a formal limit ?i — > the term can be 

neglected. In such an approximation the solution of the Schrodinger equation 
(4) is expressed by the classical flow starting from x (here < s <t) 

f = -^(Vm-.(y.) + ^A,_.(y.)) (7) 
Till o{U) terms we have 

V(t,x) = exp(^iWt{^))exp(^ - J* ^(^AW^-M + f rfi«At_,(y,))dsj ,^(yt(x)) 

If we know the trajectory (e.g., from the Hamilton equations) then we can 
compute Wt 

;y.(x, = H.(y.(x))+J^'(|(J)^ + £A,(y.)g)^ (8) 

From the correlation functions (2) of A it follows that the assumption that 

A{t, x) is x-indcpcndcnt is a good approximation for a non-rclativistic motion. 
For a wave packet of momentum p we have W = px. Hence, approximately 

f^ = -^(p+^A,_.) (9) 
ds m c 
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This is a consistent approximation because for the approximate solution of the 
Hamilton- Jacobi equation VWs(x) is space independent . Then, for A which is 
space independent we have the exact solution of eq.(5) 

As a result of the evolution in an environment of photons (which are not under 
an observation) the pure state 

ip = e-Kp{iW/h)(j} 

after an average over the states of the quantum electromagnetic field is trans- 
formed into a mixed state with the density matrix 

pt(x,x') = (?(x)^(x')) (11) 

where the average is over the electromagnetic field. 

We combine the approximation (9)-(10) of a space independent A with the 
exact Wt (8) in the following approximate expression for Wt 

^.(x) . px - _ _f_p| A.(x - ^p)dr (12) 

This expression results from a solution of the equations of motion (7) to the 
lowest (zeroth) order in e . Subsequently, Wt in eq.(8) is calculated to the first 
order in e. Inserting Wt (12) into eq.(ll) we obtain 

pt(x,x') EE exp ((ipx' - ipx)/hj(j){-K - ^p)(t)ix' - ^p)exp(-S') = 
exp f (ipx' - ipyi)/h^ ^(x - ^p)0(x' - ^p) 

exp - ^2^2j^'i /o pGth{{s - T)p/m, s - T)pdsdT (13) 
+ 2^2e2^a /o pGth(x - x' + (s - r)p/m, s ~ T)pdsdT 



'2m^c^ Jo 



/o pGth(x' - X + (s - r)p/m, s - T)pdsdT 



where exp(— S*) denotes the last factor in eq.(13). If the vacuum fluctuations 
were to be taken into account then we would need to make the replacement 
Gth — > Gf3 = Gth + *Ai?, where Af is the Feynman causal propagator (in the 
notation of Bjorken and Drell p^). 



3 The estimates of the evolution of the density 
matrix 

We perform first the integral over time in eq.(13). We denote y = x — x', 
introduce the spherical coordinates dk — dkk^dO sin 6d(j) and write S in the 
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form 

•5 - ir dkk /o" ddsineiie, k,p) (exp {pTick) - l)"' (14) 

We restrict ourselves to |p| <C mc then 
I{e, k,p) = c-2A;-2|p|2(i _ cos2 9) 

(2 (1 - cos (tck)) - 2cosky + cos (ky + ckt) + cos (ky - ckt) ^ ^^^^ 



If y is large then on the basis of eq.(15) the y-dependent terms are small as 
a function of t in comparison with other terms, because we have an additional 
oscillation in y which makes such a term negligible. The main contribution to 
5* for a large y is (we write a — cos 9) 

S - ^J^p'' J\ da{l - a^) dfcfc-i (1 - cos (tck)) (exp {(3hck) - ly' 

(16) 

For a small time t ( and a large y) we obtain 



\Pt\ 



exp ( - |;;r47P'*' C dkk (exp i(3hck) - ly' 



''dB 



(17) 



«exp (-B|l(pt/m)2 
where i? is a constant of order 1 and 

IdB = chp 

denotes the thermal de Broglie wave length at temperature T . The result (17) 
means that the decoherence (or simply the effect of the electromagnetic envi- 
ronment) is visible after a particle makes a path comparable with de Broglie 
wave length. 

Let us consider now a large i > in eq.(16). We apply the formula 

1 — cos w — w / d'-f ain{'-fw) 
Jo 

and the formula 3.911 of Gradshtein and Ryzhik [ p^ 



y du sm(au) y exp(/3M) - 1 j = — coth(— ) - — 



Then, we obtain 



S = |;™p' /o dl dk sin (tckj) (exp {f3hck) - ly' 

FP'/o^7(ikeoth(g)-i) (18) 



3 m^c^h-n 

pT ct 



^ 3 3, p^ln (^smh(^) , 
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for a large t (and a large |y|). 

Let us consider next small t together with a small y. Then, we can set 
s = T = in the argument of Gth in eq.(13). In such a case 

S ~ p^t^ r„°° dkk C ddsinO 

m-'chTr'^ JO JO (19) 

(l-cos2 6l) (l-cos(A;cos6'|y|))(exp(/37icfc)-l)"^ ^ ' 

The formula (19) is relevant only for a small y (for a large y we return to 
eq.(17)). From eq.(19) for a small y and small t we obtain 



p*(x,x') P« exp ( - ^^2£_p2t^y2 J- dhk^eM^cPk) - 1)-^) 
-exp(-5^(f)V%-i) 



(20) 



where B is a constant of order 1. Again we can conclude that the decoherence 
is visible on distances comparable to the thermal dc Broglic length. 

For a large t we neglect the quickly oscillating ^-dependent terms in eq.(15). 
Then, similarly as in the computations in eq.(18) (with a = cos9) 

S = ^j4% P^lyl !o 'i'yJ-i '^""(l - ir '^^ (°^P (/3^cfc) - 1)"^ sin(/c7Q;|y|) 
= ^|!^^p2 fida(l-a2)ln('^sinh(2^)') « fi-E^M 

(21) 

for a large y (and a large t). 

For a small |y| (and a large t ) the f-dependent terms in eq.(15) can be 
neglected. Then, the integral (14) gives 



e2 



\9\MXl\ dkk{e^pk-l)-' 
Jo 



15 Trm'^c^h 
Hence, 

|pt|«exp(-i?|^|p|2(mc)-2|y|%-j) 
At this point it is useful to recall the definition of the Wigner function 

>V(q,k) = {27rh)~^ j duexp(iku/?i)p(q + u/2, q - u/2) 

If p Ri exp(— |y^ — ipy/U) then 

W(q, k) ^ exp(--l^(p - k)2)«;(q, k) 
2an 

Such a behaviour means a localization on the classical momentum p. If (as 
in eq.(18) for a large time) pt = exp(— ipy — 5) ~ exp(— ipy — bp^t) then 
approximately 

dtPt^-b[P,[P,pt]] (22) 
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where P = —ihV is the quantum momentum operator. Then 

atWt(q,k) « -fek2>Vt(q,k) 



For N-particles with momenta Pj we consider the initial state of the form 



Then, Wt in eq.(12) is a sum of the Hamilton- Jacobi functions for each parti- 
cle. The subsequent expectation value over the electromagnetic field gives an 
exponential with a sum of pairings between the N particles. Hence, it can be 
expected that pt behaves as exp(— i?iV^t^) for a small t and as exp(— i?7V^|i|) 
for a large t with a certain constant R (note that the behaviour e'xp{—RN\t\) 
for another model has been obtained by Unruh jl^). Such a result is valid for 
N particles of equal charges under the assumption that the terms in the expo- 
nential of the form (13) have equal signs. The contributions add if the momenta 
have a distinguished direction. If the directions are random then the contribu- 
tions from different particles cancel one another. A similar cancellation takes 
place if the system is neutral , i.e., if the charges can have different signs. 
Then, in the sum in an exponential of the form (13) we shall have ejCfeP^Pfc 
multiplying Gth- Hence, there can be cancellations from different charges 
even if there is a distinguished direction of the momenta. Under an assumption 
that the contributions in the exponential do not cancel one another we obtain 
for a large distance between any pair of coordinates and a large t 



where B is a constant of order 1. Such a decay rate means that the decay is 
visible if the distance achieved by a particle (with the mean momentum p ) is 
comparable with the de Broglie wave length. The time needed for this purpose 
can be short if N is large 



4 Disappearance of the interference 

So far we have discussed the decay of pt(x, x') for varying t and y = x — x'. The 
disappearance of the off-diagonal matrix elements indicates a classical behaviour 
of quantum probabilities in the state pt- We investigate next what happens 
with the interference describing a typical quantum phenomenon. We consider a 
superposition of two wave packets 



N 






?/;(x) = exp(ip(-^^x/?i)(/)(x) + exp(ip<^^)x/?i)0(x) 
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Then, at the point x (on the screen) after an evolution through a cavity filled 
with thermal photons the probability density (|'^t(x)p) is equal to the diagonal 
part of the density matrix 

Pt = {\4't)m (23) 

For each packet we have the Hamilton- Jacobi function 

m(x) = px - - ^ / A(y,(x), s)ds (24) 
2m mc Jo 

Under the time evolution 

In our semiclassical approximation (/)(x) — >■ (/)(yt(x))). Then, for weak fields 
we neglect the dependence of the paths on the electromagnetic field, i.e., we 
consider straight lines 

and 

(2) * (2) 

m 

In this approximation the expectation value (23) is 
(<^(x-^p(2))0(x-^pW) 

exp ( - l(p(2) - p(i))x + 5S-((p(2))2 - (p(i))2)) + 

exp f i(p(2) - p(i))x - 25r((p('^)' - (p«)2))<A(x-^p(i))</.(x - ^p(2))) 
exp /o p(^)G..(^p(i) - ^p(2), . - r)p(^)dsdr 

+ /o P^^^G..(;f,p(^) - ^p(i), . - r)p(^)d.dr 

-M- /o p(^^G,.(;|p(i) - ^p(i), . - r)p(i)dsdr 

-2;^ /o P^^^G..(^p(2) - ^p(2), . - r)p(^)ci.dr) ^ p« + pf^ + p^' 

(25) 

Without detailed calculations we can obtain the behaviour for a small t. Let 
s = T = in Gt/i in eq.(25) then from eq.(3) G(?,(0,0) ~ BSji/3-^h'^c-\ It 
follows that 



exp(-i?|^/,-|(t|p(2)-p(i)|/m)2) (26) 



Hence, if the decoherence is to be visible the distance between the particles after 
time t must be of the order of the de Broglie length. The calculations for a large 
time are more involved. Let us denote 

|p(^')|^exp(-5i2) 
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We obtain for 5*12 



S,2 - 1,^\P^'^ - P^'^ P r f (exp(/3M) - (1 - cos (tck)) 
^\P'-^^ - P^'^Pct/o j;;^ dk{exp{(3nck) - l)-'sm{jckt) 

3 T^m'^c^h „ _, 



3 Trm^ 

2 'p(^)-p«P/o^d7(2i^coth(^)-^) (27) 



for a large t. Hence, the mixed term in eq.(25) is multiplied by exp(— 512) which 
decays as exp(— Such a behaviour of the probability density proves that 
the thermal photons lead to the classical addition of probabilities instead of 
the quantum addition of amplitudes showing the decoherence phenomenon in a 
physical model of an electron-photon interaction. 



5 Conclusions 

We have discussed a model of a quantum charged particle interacting with a 
quantum electromagnetic field at finite temperature. We have calculated the 
time evolution of the density matrix in a semiclassical approximation for the 
wave function and in a weak coupling approximation for the particle-photon 
interaction. In contradistinction to refs.0|^|^ we do not apply the approxi- 
mation of a linear coordinate coupling to the environment. For a small time 
and small space separations we obtain an exponential in time and space decay 
Pt ~ exp(— 6i^|x — x'p) of off-diagonal matrix elements (decoherence). For a 
large time the decay achieves its stationary value (18) and (21). The time and 
space scale is determined by the thermal de Broglie wave length. If we have a 
large number N of charged particles then the decoherence rate can increase as 
7V^ . The density matrix elements decay as pt ~ exp(— 5t) for a large time. Such 
a behaviour is in agreement with the Lindblad dynamics if the dissipative part 
of the dynamics has the form (22). Such Lindblad dynamics has been discussed 
in refs.p^p6|. The Lindblad dynamics resulting from a linear coordinate cou- 
pling to the environment (studied in |4j||^) is of a different type. It could be 
described by a replacement of the momentum operator by a position operator 
in eq.(22) (a coupling to the environment linear in the momentum as well as 
in the coordinate has been discussed by Leggett in ref.jl^). We have studied 
the interference as another typical aspect of a quantum behaviour. We have 
shown that in an environment of photons the interference disappears with an 
exponential speed (eqs.(26)-(27)). Our results suggest an arrangement for an 
experiment. Such experiments could verify the QED beyond the usual pertur- 
bative approximation as well as the principle of the wave function reduction in 
a non-selective measurement. 



9 



References 

G. R. Blumenthal and R.J. Gould, Rev.Mod.Phys.42,237(1970) 
L.S. Brown and R. S. Steinke, Am.J.Phys.65,304(1997) 
B. Dehning et al, Phys.Lett.B249,145(1990) 
A.O. Caldeira and A.J. Leggett, Physica A121,587(1983) 
W.G. Unruh and W.H.Zurek, Phys.Rev.D40,1071(1989) 
P.M.V.B. Barone and A.O.Caldeira,Phys.Rev.A43,57(1991) 
Z. Haba, Journ.Math.Phys.39,1766(1998) 
E. Jogs and H.D. Zch, Z. Phys. B59,223(1985) 
A. Stern, Y. Aharonov and Y. Imry, Phys.Rev. A41,3436(1990) 
L.H. Ford, Phys.Rev. D47,5571(1993) 

H. P. Stapp, Phys.Rev.A46,6860(1992) 

J.D. Bjorken and S.D. Drell, Relativistic Quantum Fields, McGraw- 
Hill, 1965 

W.G. Unruh, Phys.Rev.A51,992(1995) 

I. S. Gradshtcin and I.M. Ryzhik, Tables of Integrals, Sums, Series and Prod- 
ucts, Nauka, Moscow. 1971 (in Russian) 

A. Sandulcscu and H. Scutaru, Ann.Phys.(N.Y.) 173,277(1987) 
Z. Haba, Phys.Rcv.A57,4034(1998) 
A.J. Leggett, Phys.Rev.B30,1208(1984) 



10 



